Extrapolated Smoothing Descent Algorithm for Constrained Nonconvex and Nonsmooth Composite Problems*

In this paper, the authors propose a novel smoothing descent type algorithm with extrapolation for solving a class of constrained nonsmooth and nonconvex problems,where the nonconvex term is possibly nonsmooth. Their algorithm adopts the proximal gradient algorithm with extrapolation and a safe-guarding policy to minimize the smoothed objective function for better practical and theoretical performance. Moreover, the algorithm uses a easily checking rule to update the smoothing parameter to ensure that any accumulation point of the generated sequence is an (affine-scaled) Clarke stationary point of the original nonsmooth and nonconvex problem. Their experimental results indicate the effectiveness of the proposed algorithm.     

基于简化的应变梯度理论下Kirchhoff板模型边值问题的提法及其应用北大核心CSCD

考虑应变梯度和速度梯度的影响,建立薄板控制微分方程及给出其边值问题的提法,修正了前人给出的薄板角点条件.采用Levy法,给出受分布力作用下简支板的挠度及自由振动频率的解析解.通过与文献中分子动力学数据对比,验证了该文模型的有效性并提出校核材料参数的一种方法.研究结果表明,增大弹性地基和应变梯度参数可以有效提高板的等效刚度,而速度梯度参数则相反.该文提出的板的边值问题为研究薄板在复杂支撑边界及外荷载等条件提供了理论依据.同时,有望为其有限元法、有限差分法和基于能量原理的Galerkin法等数值方法提供理论依据.     

Elastic Analysis of Anisotropic Functionally Graded Rotating Disks With Non-Uniform Thicknesses北大核心CSCD

研究了任意梯度变化的变厚度各向异性转动圆盘的弹性问题.假设圆盘绕刚性轴匀速转动,其材料性能和厚度沿径向任意梯度变化.考虑圆盘在中心转轴处受位移约束,外侧自由,根据各向异性转动圆盘的平衡微分方程,得到关于径向应力的Fredholm积分方程,继而通过对Fredholm积分方程进行数值求解,得到结构的位移场和应力场.对具体梯度变化情况仅需代入相应梯度变化进行求解即可.数值算例部分,通过假设厚度、弹性模量等参数为特殊的幂函数形式,将由Fredholm积分方程求出的数值解与对应的精确解进行对比,以及针对常见的Voigt模型,将由该方法算得的数值解和ANSYS有限元计算结果进行对比,验证了该方法的准确性和精度.其次,针对Voigt模型,重点分析了厚度变化、材料性能梯度参数、各向异性度等对应力场和位移场的影响.提出了针对材料性能和厚度沿径向呈任意梯度变化的圆盘结构弹性分析方法,将为优化功能梯度圆盘的结构和材料参数、有效调整构件应力分布、提高结构安全性,提供强有力的工具;算例分析结果对功能梯度圆盘在复杂条件下的结构安全设计有重要的理论指导意义.     

修正PRP共轭梯度方法求解无约束最优化问题

基于著名的PRP共轭梯度方法,利用CG_DESCENT共轭梯度方法的结构,本文提出了一种求解大规模无约束最优化问题的修正PRP共轭梯度方法。该方法在每一步迭代中均能够产生一个充分下降的搜索方向,且独立于任何线搜索条件。在标准Wolfe线搜索条件下,证明了修正PRP共轭梯度方法的全局收敛性和线性收敛速度。数值结果展示了修正PRP方法对给定的测试问题是非常有效的。     

Fastiterative methods for least squares estimations

Least squares estimations have been used extensively in many applications, e.g. system identification and signal prediction. When the stochastic process is stationary, the least squares estimators can be found by solving a Toeplitz or near-Toeplitz matrix system depending on the knowledge of the data statistics. In this paper, we employ the preconditioned conjugate gradient method with circulant preconditioners to solve such systems. Our proposed circulant preconditioners are derived from the spectral property of the given stationary process. In the case where the spectral density functions() of the process is known, we prove that ifs() is a positive continuous function, then the spectrum of the preconditioned system will be clustered around 1 and the method converges superlinearly. However, if the statistics of the process is unknown, then we prove that with probability 1, the spectrum of the preconditioned system is still clustered around 1 provided that large data samples are taken. For finite impulse response (FIR) system identification problems, our numerical results show that annth order least squares estimator can usually be obtained inO(n logn) operations whenO(n) data samples are used. Finally, we remark that our algorithm can be modified to suit the applications of recursive least squares computations with the proper use of sliding window method arising in signal processing applications.Research supported in part by HKRGC grant no. 221600070, ONR contract no. N00014-90-J-1695 and DOE grant no. DE-FG03-87ER25037.     

Volumetric properties of ethanol–water mixtures under high temperatures and pressures

The densities of ethanol and ethanol–water mixtures were measured with a vibrating tube densimeter at 25.0, 50.0 and 75.0 °C in the pressure range from 0.10 to 40.00 MPa. Densities were correlated using an empirical model. Partial molar volumes, excess molar volumes, isothermal compressibilities, cubic expansion coefficients and internal pressures were calculated from obtained densities. This study reports the dependence of densities, partial molar volumes, excess molar volumes, isothermal compressibilities, cubic expansion coefficients and internal pressures on composition, temperature and pressure.     

High-performance liquid chromatography with continuous gradient elution in screening for unknown metabolites of butoprozine directly in dog bile

This paper describes the screening for metabolites of butoprozine, a new anti-anginal drug, in dog bile by means of reverse-phase HPLC. Although it does involve a simple clean-up step to remove a substantial amount of endogenous bile compounds, this screening method nevertheless avoids extraction of metabolites and thus allows all metabolites to be introduced into the chromatographic system. A single run of 100 minutes from 100% water to 100% methanol in a linear gradient effects adequate separation of the great majority of metabolites without interference from remaining endogenous compounds. Two methods of differentiating between metabolite peaks and endogenous peaks have been worked out. The first one makes use of ¹⁴C-labeled butoprozine by measuring the amount of radioactivity in the column effluent while simultaneously recording the UV absorbance. The second method compares continuous gradient chromatograms of bile recorded before and after butoprozine administration under very similar conditions. The latter method can be applied to both radioactive and non-radioactive materials.     

Exploring high pressure structural transformations,electronic properties and superconducting properties of MH2 (M = Nb,Ta)

The high-pressure structures and properties of MH₂ (M = Nb, Ta) are explored through an ab initio evolutionary algorithm for crystal structure prediction and first-principles calculations. It is found that NbH₂ undergoes a phase transition from a cubic Fm$\overline{3}$m structure with regular NbH₈ cubes to an orthorhombic Pnma structure with fascinating distorted NbH₉ tetrakaidecahedrons at 48.8 GPa, while the phase transition pressure of TaH₂ from a hexagonal P6₃mc phase with slightly distorted TaH₇ decahedron to an orthorhombic Pnma phase with attractive distorted TaH₉ tetrakaidecahedrons is about 90.0 GPa. Besides, the calculated electronic band structure and density of states demonstrate that all of these structures are metallic. The Poisson’s ratio, electron localization function, and Bader charge analysis suggest that these phases possess dominant ionic bonding character with the effective charges transferring from the metal atom to H. From our electron–phonon calculations, the calculated superconducting critical temperature T_c of the Pnma-NbH₂ is 6.903 K at 50 GPa. Finally, via the quasi-harmonic approximation method, the phase diagrams at pressure up to 300 GPa and temperature up to 1000 K of MH₂ (M = Nb, Ta) are established, where the transition pressure of Fm$\overline{3}$m-NbH₂ → Pnma-NbH₂ and P6₃mc-TaH₂ → Pnma-TaH₂ were found to decrease with increasing temperature.     

Numerical examination of performance of some exchange-correlation functionals for molecules containing heavy elements

The performance of 17 exchange-correlation functionals for molecules containing heavy elements are numerically examined through four-component relativistic density DFT calculations. The examined functionals show the similar accuracy as they do for the molecules containing light elements only except for bond lengths. LDA and OP86 produce good results for bond lengths and frequencies but bad bond energies. Different functionals do not show much different performance for bond energies except LDA. BP86 and GP86 produce results with average accuracy while LYP does not perform well. Although encouraging results are obtained with functional B97GGA-1, other heavily parameterized and meta-GGA functionals do not produce impressive results.     

Nondifferential optimization via adaptive smoothing

The problem of minimizing a nondifferential functionx f(x) (subject, possibly, to nondifferential constraints) is considered. Conventional algorithms are employed for minimizing a differential approximationf off (subject to differentiable approximations ofg). The parameter is adaptively reduced in such a way as to ensure convergence to points satisfying necessary conditions of optimality for the original problem.This research was supported by the UK Science and Engineering Research Council, the National Science Foundation under Grant No. ECS-8121149, and the Joint Services Electronics Program, Contract No. F49620-79-C-0178.